Transport criticality in triangular lattice Hubbard model 
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We study electric transport near the Mott metal-insulator transition. Optical conductivity of the 
half-filled Hubbard model on a triangular lattice is calculated based on a cellular dynamical mean 
field theory including vertex corrections inside the cluster. By investigating the spectrum at low 
frequencies, we find that a Drude peak on the metallic side smoothly connects to an "ingap" peak 
on the insulating side. The optical weight of these peaks exhibits a critical behavior with power-law 
near the Mott critical end point, \D — D*\ oc \U — [7*| 1//<s . We find that the critical exponent 1/5 
differs from the exponents in the thermodynamics. 

PACS numbers: 71.27.+a, 71.30.+h, 72.10.-d 



Mott transition is one of the central topics in strongly 
correlated electronic systems. This occurs when Coulomb 
repulsion predominates over electric kinetic energy. Mott 
transition takes place inside the paramagnetic region and 
magnetic frustration plays an important role to real- 
ize it. A classic three-dimensional example is Cr-dopcd 
V203fi and quasi- two-dimensional K-type organic com- 
pounds have been intensively studied recently— Their 
phase diagram in the parameter space of temperature 
T and (chemical) pressure P has insulating and metallic 
phases which are separated by a line of first-order Mott 
transition, and this phase boundary terminates at an end 
point. This is similar to the liquid-gas transition in clas- 
sical liquids^ and critical behaviors are expected for vari- 
ous properties around the Mott critical end point. In re- 
cent works, critical behaviors in thermodynamic quanti- 
ties have been confirmed near the Mott critical end point 
and the criticality has been discussed actively^ How- 
ever, studies on nonequilibrium or transport properties 
are limited, although the singularity is most prominent 
in electric conductivity^ Putative criticality in electric 
transport is the main issue of this paper. We will demon- 
strate its existence at the Mott critical point and analyze 
its singularity by numerical approach. 

Recently, two experimental studies have reported a 
critical behavior of dc-conductivity oq. Limclettc et 
al. performed a scaling analysis for <jq of (Vi_ x Cr x )203 
near the Mott critical end point, which indicates 
the three-dimensional Ising universal class (UC). 6 
However, using the similar analysis, Kagawa et 
al. found that the critical behaviors in the quasi-two- 
dimensional anisotropic triangular lattice compound k- 
(ET)2Cu[N(CN)2]C1 are not classified as any conven- 
tional UC'sii There are several proposals to explain the 
latter. Imada et al. suggested that the unconventional 
UC is due to a marginal quantum critical point in the 
two-dimensional system^ Another theory by Papaniko- 
laou et al. proposed that the different exponents come 
from the fact that conductivity corresponds to energy 
density of the Ising model in addiition to magnetization. 9 

Optical conductivity, <t(w), also shows peculiar behav- 
iors in the K-type organic family, depending on geomet- 




FIG. 1: (Color online) Optical conductivity <r(w) for various 
U at T = 0.10 near the critical point. Inset is the U-T phase 
diagram. Diamond (circle) represents the boundaries of the 
coexisting phase (crossover region). The line of the first order 
transition is a guide for eyes. 

rical frustration in each material. Strongly frustrated 
Mott insulator k-(ET) 2 Cu2(CN) 3 exhibits a spin liquid 
behavior, and its optical conductivity does not show a 
clear gap but a smooth decay toward zero frequency^ 
In contrast, k-(ET) 2 Cu[N(CN) 2 ]C1 is less frustrated due 
to a distorted triangular structure, and its <j{uj) exhibits 
a clear gap. 11 ' 12 

In this paper, we numerically study the "Mott criti- 
cality" in transport properties of a strongly correlated 
electronic system on a geometrically frustrated lattice. 
In particular, we focus on optical conductivity near the 
Mott critical end point and examine its behaviors. 

The model to study is the one-band Hubbard Hamil- 
tonian on an isotropic triangular lattice at half filling, 

H = -t 4a c Ja + U J2 ni t n 4 - M C^Cur, (1) 

(i,j},tr i i,a 

where t is the nearest-neighbor hopping amplitude, U 
is the on-site Coulomb repulsion and /j, is the chemical 
potential to tune electron density. ci a is an electron an- 
nihilation operator at site i with spin a and rii a — c irj Cic. 
To calculate optical conductivity, we need single- and 
two-electron Green's functions. We compute them by 
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the cellular dynamical mean field theory (CDMFT)^ to 
take into account both strong electronic correlations and 
geometrical frustration. We map Hamiltonian (p} onto a 
three-site cluster model coupled to an effective medium 
determined self-consistcntly. The cluster Green's func- 
tions are computed by using the continuous-time Quan- 
tum Monte Carlo (CTQMC) method based on the strong 
coupling expansion^ 

Optical conductivity <t{oj) is calculated from the 
current-current correlation function Xqi 1 ^) with fre- 
quency oj and wave vector q — ► based on the Kubo 
formula. In the single-site DMFT approach, Xq( w ) 
is usually calculated by convoluting the single-electron 
Green's function^ To take into account correlation ef- 
fects, we proceed further beyond the standard formu- 
lation and include vertex corrections^ inside the clus- 
ter, which is a big challenge in numerical computations. 
Our study with CDMFT is one of the first achievements 
for conductivity and the only preceding work studied a 
square lattice system and employed the dynamical cluster 
approximation^ We will explain some details later and 
first report results. In what follows, we normalize u{uj) 
by the unit of e 2 /fi, where e is the elementary charge and 
H is the reduced Planck's constant and U and tempera- 
ture T are in units of t. 

We start with summarizing the variation of ct(w) near 
the Mott critical end point. The inset of Fig. [T] shows the 
phase diagram determined by our calculation of double 
occupancy in the parameter space of U and TJ£ A line of 
first-order Mott transition separates the metallic phase 
from the insulating phase, and this line terminates at 
the end point, U* ~ 9.4 and T* ~ 0.1. The main panel 
shows the variation of <t(uj) with increasing U at T — 0.10 
fixed, which is nearly T* , and we have checked there is 
no hysteresis. For U < 9.4, cr(u>) shows a Drude peak at 
low w, indication of the metallic state, and an incoherent 
broad peak around oj ~ U. For U > 9.5, the Drude peak 
disappears quickly and its absence is a characteristic of 
insulator. We also calculated <t(lu) at various T's. At 
the lower T = 0.09, o(uj) shows a jump and hysteresis, 
corresponding to the first-order transition. At the higher 
T = 0.15, cr(ui) shows a smooth crossover from metal 
to insulator. We will investigate expected singularity of 
c(uj), and analyze its dependence on U at T = 0.10 fixed. 

We first discuss the metallic side. We fit the low energy 
peak by a simple Drude formula <r{uS) — Re[Do/(— iu; + 
1/t)] and analyze the [/-dependence of Do and 1/r. Dq 
is a Drude weight, which is inversely proportional to the 
effective mass to*, and 1/r is the transport scattering 
rate. One may expect that 1/r diverges with approaching 
U*. Instead, we find that 1/r is almost constant and D 
drops drastically near the Mott transition point, implying 
the increase of to* as shown in Fig. [5] (a). These results 
are consistent with an extended Drude analysis including 
w-dependence in m* and 1/r. 

On the insulating side, in addition to a broad peak 
around uj ~ U , cr(u>) has a small peak at a lower ui. The 
former persists from the metallic side and comes from the 




2 4 6 



CO 

FIG. 2: (Color online) (a) Dependence on U of the weights 
of Drude peak Do, an ingap peak Dig, and total weight 
D = Do + Dig- Inset is [/-dependence of transport scat- 
tering rate 1/r on the metallic side, (b) Ingap peak <tig(w) 
on the insulating side for U > 9.4. An ingap peak exists at 
cj ~ 1 and a small Drude peak coexists at U = 9.4. 

excitations to the Hubbard band. Analysis reveals that 
its peak position shifts from the strong coupling limit 
u> ~ U by a finite amount ~ —3, and this is due to 
the motion of doublon-holon pairs. This Hubbard "gap" 
in <r(w) is not so clear and the incoherent peak is well 
fitted by a simple Gaussian form, particularly for to below 
the peak position. We subtract the Hubbard part from 
the total cr(u>) and extract the low-w part, ctig(w). The 
results are shown in Fig. [2] (b). They show a peak below 
u) ~ 2, and with approaching U*, the peak position shifts 
toward lower u> and the intensity grows quite drastically. 
We will call this structure an ingap (IG) peak. The IG 
peak comes from low-energy peaks with small intensity 
in the density of states near the Mott transition point. 

Let us start a detailed analysis of ct(uj) near U* and 
examine whether it exhibits a singularity. Since er(cj) is 
not derived from free energy, protocol of its analysis is 
not standardized like thermodynamic criticality and we 
try several quantities as a scaling variable. The simplest 
choice is dc-conductivity <7q — er(0) and it was used in 
experimental studies^™ However, it turns out that this is 
not convenient on the insulating side, as we will explain 
later, do is also very sensitive to impurity scatterings, 
which perturb intrinsic behavior. An alternative choice 
of scaling variable is the Drude weight Do on the metallic 
side and this is more robust against impurity scatterings. 
On the insulating side, however, Do = 0, but there exists 
an IG peak at low-w, which evolves into a Drude peak 
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FIG. 3: (Color online) Log-log plot for the scaling analysis of 
the optical weight D. Scaling functions determined by fitting 
for 4 data points on each side are shown by straight lines. 
Inset shows the variation of the exponents with relaxing the 
fitting error by the amount specified. 

on the metallic side. Thus, it is natural to use its weight 
D\g = (2/7r) L cIlj<jig(ll>) as a counterpart of Dq. 

Figure[3](a) shows the [/-dependence of D = D + Diq- 
D changes continuously with U but exhibits a singularity; 
a divergent slope near U = 9.4. Following the idea of 
thermodynamic criticality, we fit the curve with a power- 
law scaling function, 

\D-D*\ = A±\U -U*\ 1/s± , (2) 

where the index — (+) denotes the metallic (insulating) 
side. Figure [3] shows the result of scaling analysis using 4 
points on each side around the critical point. The fitting 
is successful on both sides and D does exhibit the critical 
behavior. The critical point is determined simultaneously 
in this fitting (U*,D*) = (9.3999±10~ 4 , 0.11±0.01). The 
error is estimated by relaxing the error minimization in 
the fitting by 5 %. 

The most important result is the critical exponents, 
(l/<J_,l/5+) = (0.15 ±0.02, 0.13 ±0.02). While 1/6+ = 
l/<5_ as expected within the accuracy, this value is un- 
expected. See the inset in Fig. [3] In the analysis of ex- 
perimental data^ii the fundamental assumption is that 
deviation in pressure AP from the critical point and dc- 
conductivity A<7o should be interpreted as magnetic field 
h and magnetization m in the Ising ferromagnct. There- 
fore, the exponent obtained from conductivity is believed 
to coincide with the one relating m oc Ik] 1 / 6 at the critical 
temperature in the Ising system. The value 1 /6 (Ising) is 
1/15 and 1/4.8 for the two and three-dimensional cases, 
and 1/3 in the mean field theory, which becomes exact in 
four and higher dimensions^ The exponent determined 
in our analysis does not agree with any of these values. 
Our result is unexpected and interesting, since critical- 
ity in thermodynamics is well described by the Ising UC. 
In the Mott transition, double occupancy (n^n^) plays 
the role of order parameter, and the DMFT calculation 
shows that its dependence on U and T is well fitted by the 
scaling function with the mean filed exponents ft = 1/2 



and 1/5 = 1/3. Our result is also different from the one 
by Papanikolaou et al£ Therefore, the present result im- 
plies the breakdown of the analogy with thermodynamic 
criticality and we need a new description for transport 
properties beyond the one based on free energy, since 
they are nonequilibrium phenomena. 

One needs some care in the scaling analysis. We 
first tried the above scaling analysis with data points 
a little more apart from the critical point and the ob- 
tained exponents were distinct between the two sides, 
1/6- = 0.27 ± 0.03 and 1/6+ = 0.11 ± 0.03. We 
get to obtain 5—~6+ after a few more data points are 
added closer to the critical point. This indicates that 
the true critical region is rather narrow and subleading 
terms are not negligible outside it. The critical region is 
—0.12 < U — U* < 0.11, which is determined such that 
the subleading term contributes less than 2 % there. 

It is also important to note that we obtain the same 
exponent on the insulating side only when D is chosen 
as a scaling variable. We performed the same analysis 
for dc-conductivity <jq. Its suppression on the insulating 
side is more prominent than for D, and the obtained 
exponents are clearly distinct, (l/6-,l/6+) = (0.16 ± 
0.04,0.07 ± 0.02), but the value on the metallic side is 
close to the one determined from D. This robustness 
is also consistent with the fact that another scaling with 
using only Do, Da = A(U* — U) 1 ^ , also leads to a similar- 
value 1/6 = 0.13±0.03. Scaling analysis of experimental 
data was performed only for the metallic side^ii and it 
is useful to check by experiments whether the insulating 
side has the same exponent. 

Finally, let us briefly explain our algorithm of including 
the vertex corrections^ in CDMFT. Based on the Kubo 
formula, a(w) is defined as u(uj) = — ilim q ^o[Xq( w ) ~ 
Xq(0)]/w. Here, the current-current correlation function 
is described in Matsubara space as, 

J I ■ \ 2 O.cr/. \ 

Xq = o( lW «) = W kX k 

+ VkVk'Xk CT (^™) r kk'(* w n)Xk' <T (* w «),(3) 

where X^i^n) = -jjG^(iei)G^,(iei + iuj n ) and G^(iei) 
is the Green's function of electron with wave vector 
k and spin a. The above equation is exact except 
that the dependence on a and e;/ has been averaged 
over in the full vertex r££, . We use the convention 
of summing over repeated indices. N is total num- 
ber of sites and is the rr-component of velocity 

Vk = 2t(sinfc x + sinlf cos-^!^). Spin susceptibility has 
been calculated including vertex corrections already in 
single-site DMFT approaches,— while vertex corrections 
in current correlations can be included only in the cluster 
generalization. To obtain the full vertex r££, (iu) n ), we 
first calculate two-electron Green's functions in the 
cluster x^' 7 ,(r) = (cL(r)c^(r)4 CT ,(0)c^(0)) by 
CTQMC, where a-6 denote sites in the cluster. From 
this, we obtain the irreducible vertex in the cluster 
by solving the Bethe-Salpeter equation Xap-ysi^n) = 

Xapysfan) + XaJa'^( iw ") J a'X'5'( iw ™)X7'X«( ia; ™)' 
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FIG. 4: (Color online) Effects of vertex corrections on optical 
conductivity o~{uj). 

where Xda^si^n) denotes the free part calculated 
from G^(iei). We then calculate the lattice irreducible 
vertex I%£,{iu n ) = I^ T4 (*w„)e*( r «- r ^^< r -T- r »>, 
and obtain T^,{iuj n ) via the Bethe-Salpeter equation, 

(tw„) = {iu> n ) + [iu n ) X ^" (iun)T&£ (tw„). 
We finally calculate u{ui) by means of analytic contin- 
uation iu) n — > 10 + iO by using the maximum entropy 
algorithm^ 

Figure [4] presents the effects of vertex correc- 
tions; crb u b(w) is the results without vertex corrections 
obtained from the first term in Eq. ([3]), and the contribu- 
tion of the vertex corrections is <j v (lu) = <r(u>) — o"b u b(w). 
The data for U = 8.5 are typical result on the metal- 
lic side. The vertex correction changes the Drude peak 
into more coherent; a larger weight and a narrower 
width. On the insulating side U = 9.5, an IG peak 
is enhanced and the Hubbard peak slightly shifts to- 
ward a higher u> by the vertex correction. On both 
sides, the total spectrum shifts toward a lower to. For 
the scaling analysis of D, the results determined by 
cr bub (u;) are (U*,D*) = (9.3997 ± 10" 4 , 0.09 ± 0.01) and 



1/5+) = (0.16 ±0.03, 0.12 ±0.02). Detailed anal- 
ysis shows that subleading terms have larger contribu- 
tion compared with the case including vertex corrections. 
Consequently, the critical regions are narrower, particu- 
larly on the insulating side: -0.10 < U-U* < 0.02. This 
results in larger errors in the obtained value of the expo- 
nent l/5±. Introducing local vertex correction partially 
recovers the otherwise violated charge conservation in 
DMFT, and enhances coherence in charge transport. Our 
results suggest that the vertex corrections also suppress 
subdominant modes in the critical behavior of transport 
and the criticality becomes more apparent. 

In this paper, we have studied optical conductivity 
cr(ui) near the Mott critical end point in the triangular- 
lattice Hubbard model by CDMFT combined with a 
new algorithm for including vertex corrections. We have 
found that o~(oS) has an "ingap" peak within the Hub- 
bard gap on the insulating side and that it continuously 
evolves into Drude peak on the metallic side. We ex- 
amined the possibility of criticality in a(uj) and demon- 
strated that the weight of Drude and ingap peak exhibits 
a power-law singularity at the Mott critical point. Al- 
though our calculations are a dynamical mean-field ap- 
proximation, its critical exponent differs from the pu- 
tative mean-field exponent and also from that of Ising 
universality class in any dimension. This indicates a 
nonequilibrium nature of transport criticality. 

The authors are grateful to T. Ohashi and H. Kusunose 
for helpful discussions. The present work is supported by 
MEXT Grant-in- Aid for Scientific Research on Priority 
Areas "Novel States of Matter Induced by Frustration" 
(No. 19052003), and by Next Generation Supercomput- 
ing Project, Nanoscience Program, MEXT, Japan. Nu- 
merical computation was performed with facilities at Su- 
percomputer Center in Institute for Solid State Physics 
and Information Technology Center, University of Tokyo. 



1 D. B. McWhan, et al, Phys. Rev. Lett. 27, 941 (1971). 

2 K. Kanoda, J. Phys. Soc. Jpn. 75, 051007 (2006). 

3 H. E. Stanley, Introduction to phase transition and critical 
phenomena (Oxford University Press, 1971). 

4 G. Kotliar, E. Lange, and M. J. Rozenberg, Phys. Rev. 
Lett. 84, 5180 (2000). 

6 M. Imada, Phys. Rev. B 72, 075113 (2005). 

6 P. Limelette et al, Science 302, 89 (2003). 

7 F. Kagawa, K. Miyagawa, and K. Kanoda, Nature 436, 
534 (2005). 

8 M. Imada, Phys. Rev. B 72, 075113 (2005). 

9 S. Papanikolaou et al, Phys. Rev. Lett. 100, 026408 
(2008). 

10 I. Kezsmarki et al, Phys. Rev. B 74, 201 101 (R) (2006). 

11 K. Kornelsen et al, Solid State Commun. 81, 343 (1992). 

12 T. Sasaki et al, Phys. Rev. B 69, 064508 (2004). 

13 M. J. Rozenberg et al., Phys. Rev. Lett. 75, 105 (1995). 



14 P. Nozieres, Theory of Interacting Fermi Systems (Addison 
Wesley, 1964). 

15 Nan Lin, Emanuel Gull, and A. J. Millis, Phys. Rev. B 80, 
161105(R) (2009). 

16 Y. Furukawa, Master's thesis, Kyoto University (2010). 

17 S. V. Dordevic et al, Phys. Rev. Lett. 86, 684 (2001). 

18 L. Degiorgi, F. B. B. Anders, and G. Griiner, Eur. Phys. 
J. B 19, 167 (2001). 

19 G. Kotliar, S. Y. Savrasov, G. Palsson, and G. Biroli, Phys. 
Rev. Lett. 87, 186401 (2001). 

20 P. Werner et al, Phys. Rev. Lett. 97, 076405 (2006). 

21 V. L. Pokrovsky and A. L. Talapov, Phys. Rev. Lett. 42, 
65 (1979). 

22 M. Jarrell, Phys. Rev. Lett. 69, 168 (1992). 

23 M. Jarrel and J.E. Gubernatis, Phys. Rep. 269, 133 (1996). 



